Preprint typeset in JHEP style - HYPER VERSION 



arXiv:1006.0470 [hep-ph] 
IFIC/10-06 
FTUV/10-0602 



Charged- Higgs phenomenology in the Aligned 
two-Higgs-doublet model 



Martin Jung, Antonio Pich, Paula Tuzon 

Instituto de Fisica Corpuscular (IFIC),CSIC-Universitat de Valencia 
Apartado de Correos 22085, E-46071 Valencia, Spain 



E-Mail: jung@ific.uv.es 
E-Mail: pich@ific.uv.es 



E-Mail: tuzon@ific.uv.es 



Abstract: The alignment in flavour space of the Yukawa matrices of a general two-Higgs- 
doublet model results in the absence of tree-level flavour-changing neutral currents. In 
addition to the usual fermion masses and mixings, the aligned Yukawa structure only 
contains three complex parameters which are potential new sources of CP violation 
Q . For particular values of these three parameters all known specific implementations of 
the model based on discrete Z2 symmetries are recovered. One of the most distinctive 
features of the two-Higgs-doublet model is the presence of a charged scalar . In this 
work, we discuss its main phenomenological consequences in flavour-changing processes at 
low energies and derive the corresponding constraints on the parameters of the aligned 
two-Higgs-doublet model. 
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1. Introduction 

The simplicity of the idea of including one additional Higgs doublet to the Standard Model 
(SM) and the versatility of the resulting phenomenology are the main ingredients that 
have made the two-Higgs-doublet model (2HDM, see e.g. [Q, [|] and references therein) so 
interesting. In the most general version of the model, the fermionic couplings of the neu- 
tral scalars are non-diagonal in flavour and, therefore, generate unwanted flavour-changing 
neutral-current (FCNC) phenomena. Different ways to suppress FCNCs have been devel- 
oped, giving rise to a variety of specific implementations of the 2HDM. The simplest and 
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most common approach is to impose a Zi symmetry forbidding all non-diagonal terms in 
the Lagrangian jl. Depending on the charge assignments under this symmetry, the model 
is called type I g g, II |, % X and Y g |, [R], |IJ ||, ||] or inert p, ||, 0, ||. In 
these types of models with natural flavour conservation the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix |0|, [2(| is the only possible source of CP violation. An- 
other possibility is to assume particular Yukawa textures which force the non-diagonal 
Yukawa couplings to be proportional to the geometric mean of the two fermion masses, 
gij cx ^/m~m~ the so-called type III 2HDM (H, ||, ||, |4|. 

Our work focuses on the recent suggestion []l]] to enforce the alignment in flavour space 
of the Yukawa couplings of the two scalar doublets, which guarantees the absence of tree- 
level FCNC interactions. The Yukawa structure of the resulting aligned two-Higgs-doublet 
model (A2HDM) is fully characterized by the fermion masses, the CKM quark mixing 
matrix and three complex parameters (/ = u,d,l), whose phases are potential new 
sources of CP violation The usual models based on Z2 symmetries are recovered for 
particular (real) values of these three parameters. The A2HDM provides a more general 
setting to discuss the phenomenology of 2HDMs without tree-level FCNCs, leaving open 
the possibility of having additional CP-violating phases in the Yukawa sector beyond the 
CKM-matrix one. 

The presence of a charged scalar is one of the most distinctive features of an 
extended scalar sector. In the following we analyze its phenomenological impact in low- 
energy flavour-changing processes within the A2HDM, and constrain the three complex 
parameters <?/ with present data on different leptonic, semileptonic and hadronic decays. 
We proceed as follows: the formulation of the general 2HDM is recalled in section |2|, 
where the aligned condition is implemented and the resulting Yukawa structure discussed. 
Section |3| explains our statistical treatment of theoretical uncertainties and compiles the 
inputs used in our analysis. The phenomenological consequences of having a charged scalar 
field are analyzed next, process by process, extracting the corresponding constraints on the 
new-physics parameters In section || we discuss the constraints derived from tree- level 
leptonic and semileptonic decays, while section |5] describes the information obtained from 
loop-induced processes. Finally, we give our conclusions in section ||. Some technical 
aspects related to AF = 2 transitions have been relegated to the appendix. 



2. Aligned two-Higgs-doublet model 

The 2HDM extends the SM with a second Higgs doublet of hypercharge Y = \. The 
neutral components of the scalar doublets (f> a (x) (a = 1,2) acquire vacuum expectation 
values (VEVs) that are, in general, complex: (0|<^(a;)|0) = (0, v a e l9a ). Through an 
appropriate U(l)y transformation we can enforce Q\ = 0, since only the relative phase 
9 = 62 — 61 is observable. The combination v = \J v\ + v\ ~ (y/2 Gf)^ 1 ^ 2 = 246 GeV plays 
the role of the SM VEV when generating the gauge boson masses. 

A global SU(2) transformation in the scalar space (cpi,cp2) takes us to the so-called 
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Higgs basis (<E>i,$2)> where only one doublet acquires a VEV: 




v\ v 2 

V 2 —Vl 



(2.1) 



In this basis, the two doublets are parametrized as 



$1 



G+ 

-L {v + Sl + iG°) 



c[> 9 



H+ 
^(S 2 + *S 3 



(2.2) 



where G ± and G° denote the Goldstone fields and (H+) = (G+) = (G°) = (Si) = 0. 
The five physical scalars are given by two charged fields H^(x) and three neutral ones 
<Pi(x) = {h(x), H(x), A(x)}, which are related to the Si fields through an orthogonal 
transformation y^(x) = IZijSj(x). The form of 7£y depends on the scalar potential, which 
could violate CP in its most general version; in that case the resulting mass eigenstates do 
not have a definite CP parity. 

The most general Yukawa Lagrangian of the 2HDM is given by 

c Y = - {QUri^i + r 2 <h) 4 + O' L (Ai0i + a 2 2 ) u ' r + Z^n^i + n 2 2 ) 4} + h.c. , 

_ (2.3) 

where Q' L and L' L are the left-handed quark and lepton doublets, respectively, and (f> a (x) = 
iT 2 (p* a (x) the charge-conjugated scalar doublets with Y = —5. All fermionic fields are 
written as A^-dimensional vectors and the couplings T a , A a and Il a are Nq x Nq complex 
matrices in flavour space, Nq being the number of fermion generations. Moving to the 
Higgs basis, the Lagrangian reads 



C 



V2 



Y 



with 



v 

-f h.c. , 



{Q' L {M' d <5>i + Yfi 2 ) d' R + Q' L (M'J&i + Y$ 2 ) u' R + Li(M/*i + Y{$ 2 ) I'r] + 



M^ = — (v 1 r 1 + v 2 r 2 e* 



M' u = -±= (v 1 A 1 + v 2 A 2 e- i< -' 



Ml 



V2 



villi + 



V2 



v 2 Fi 
- v 2 Ai 
uin 2 e ie - v 2 Il 



(2.4) 

(2.5) 
(2.6) 
(2.7) 



In general, the complex matrices Mj and Y* (/ = d, u, I) cannot be simultaneously diago- 
nalized. Thus, in the fermion mass-eigenstate basis, with diagonal mass matrices Mf, the 
Yukawa-coupling matrices Yj remain non-diagonal giving rise to FCNC interactions. 

The unwanted non-diagonal neutral couplings can be eliminated requiring the align- 
ment in flavour space of the Yukawa matrices [0. It is convenient to implement this 
condition in the form: 



A 2 = C u e i6 A 1 , 



n. 



6e-^ni 



(2.8) 
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Type I 


(oo, oo, oo) 


cot 8 


cot/3 


cot/3 


Type II 


(0,00,0) 


— tan/3 


cot/3 


— tan/3 


TypeX 


(oo, oo, 0) 


cot/3 


cot/3 


— tan /3 


Type Y 


(0, oo, oo) 


— tan Q 


cot/3 


cot/3 


Inert 


(tan B, tan B, tan B) 












Table 1: Limits on £f that recover the different Z2 models and corresponding values. 



where are arbitrary complex parameters. The proportionality of the matrices Yj and 
M'i guarantees that all FCNC couplings vanish at tree level: 

c . tan 8 

Y dtl = ^,iM dtl , Y U = ^M U , c/ = — ^ tan/3 . (2.9) 

! + £_/• tan p 

In the A2HDM the mass-eigenstate Yukawa Lagrangian reads ||] 

C Y = -^H + (x) {u(x) [ ft FM d ^ - ^ M+VTzJ + q P(x)M z P R /(x)} - 

- " E f/ ? ^ /» MfVnfix) + h.c. , (2.10) 

where V denotes the CKM matrix, "P^z, = 1 ^ 75 are the right-handed and left-handed 
projectors and the couplings of the neutral scalar fields are given by: 



v% = n * + ( ni2 + * n ^ ^ > y$ = n * + (^2 - i n i3 ) < . (2.11) 

Some conclusions can be drawn from (T2.10|) . In the A2HDM all fermionic couplings 
to scalars are proportional to the corresponding fermion masses and the neutral-current 
interactions are diagonal in flavour. The only source of flavour-changing interactions is 
the CKM matrix in the quark charged current, while all leptonic couplings are diagonal 
in flavour because of the absence of right-handed neutrinos in our framework, which could 
however easily be included. There are only three new parameters c/, which encode all 
possible freedom allowed by the alignment conditions; these couplings satisfy universality 
among the different generations, i.e. all fermions with a given electric charge have the same 
universal coupling c/. The three parameters are also invariant under global SU(2) trans- 



formations of the scalar fields 4> a — > <p' a = U a b4>i, |ff5 |, i.e. they are scalar-basis independent. 
Taking the particular values shown in table 1, the different models based on Z2 symmetries 
are recovered, with a single scalar doublet coupling to each type of right-handed fermions 
Q. Finally, it should be pointed out again that c/ are arbitrary complex numbers, opening 
the possibility of having new sources of CP violation without tree- level FCNCs. 

2.1 Quantum corrections 

Quantum corrections induce some misalignment of the Yukawa coupling matrices, generat- 
ing small FCNC effects suppressed by the corresponding loop factors. However, the special 
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structure of the A2HDM strongly constrains the possible FCNC interactions [||. Obvi- 
ously, the alignment condition remains stable under renormalization when it is protected 
by a Z2 symmetry f26|| , i.e. for the particular cases indicated in table |]. In the most gen- 
eral case loop corrections do generate some FCNC effects, but the resulting structures are 
enforced to satisfy the flavour symmetries of the model. The Lagrangian of the A2HDM is 
invariant under flavour-dependent phase transformations of the fermion mass eigenstates 
(/ = d,u,l,u, X = L,R, a%' = d- L ): 

fxix) e ia '' X P x {x) , Va -> e"';"l;, e , M f:ij -> e^M^ e . (2.12) 

Owing to this symmetry, lepton-flavour- violating neutral couplings are identically zero to all 
orders in perturbation theory, while in the quark sector the CKM mixing matrix remains the 
only possible source of flavour-changing transitions. The only allowed local FCNC struc- 
tures are of the type u L V{M d M ] d ) n V\M u Mt) m M u u R , d L V\M u Mt) n V{M d M j d ) m M d d R , 
or similar structures with additional factors of V, W and quark mass matrices |l]]. There- 
fore, at the quantum level the A2HDM provides an explicit implementation of the popular 



Minimal Flavour Violation scenarios M, M, pfj, 31, p2j, but allowing at the same time 



for new CP-violating phases. 1 Structures of this type have been recently discussed in fl33| 



Using the renormalization- group equations [26, 34 1, one can easily check that the one 



loop gauge corrections preserve the alignment while the only FCNC structures induced by 
the scalar contributions take the form [35|: 

£f cnc = 4^3 ( 1 + Vd) x 

x ]T <p!(x) { {K, a + i K i3 ) - ?J [d L V*M u Ml VM d d R ] - (2.13) 

i 

-(Ka-ina)kd-&) [u L VM d M d V^ M u u R ]} + h.c. 

As they should, these FCNC effects vanish identically when q = <z u (Z2 models of type 
I, X and Inert) or q = —1/?^ (types II and Y). The leptonic coupling q does not induce 
any FCNC interaction, independently of its value; the usually adopted Z2 symmetries are 
unnecessary in the lepton sector. Assuming the alignment to be exact at some scale //o, 
i.e. C(/io) = 0, a non-zero value for the FCNC coupling, C(fi) = — log (/////q), is generated 
when running to a different scale. 

The numerical effect of these contributions is, in any case, suppressed by m q m 2 q ,/v^ and 
quark-mixing factors. This implies an interesting hierarchy of FCNC effects, avoiding the 
stringent experimental constraints for light-quark systems, while allowing at the same time 
for potential interesting signals in heavy-quark transitions. Obviously, the most relevant 
terms in ( 2.13|) are the sib R and cit R operators. The sib R term induces a calculable 



contribution to B®~B® mixing through tp® exchanges, which modifies the mixing phase 
and could explain the like-sign dimuon charge asymmetry recently observed by DO p6|| . 
Tree-level scalar exchanges from FCNC vertices have been already suggested as a possible 



1 Minimal flavour violation within the context of the Type II 2HDM has been discussed in |^|. This 
reference didn't consider the possibility of incorporating new CP-violating phases. 
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explanation of the DO measurement [[37J. We defer the phenomenological analysis of the 
FCNC operator ( |2.13 ) to a future publication |3j|, where the neutral sector of the A2HDM 



will be studied in detail. In the present paper we will concentrate in the phenomenology 
of the charged-scalar Yukawa Lagrangian ( [2.10 ). 



3. Inputs and statistical treatment 

In the following sections we will analyze the most important flavour-changing processes that 
are sensitive to charged-scalar exchange and will try to constrain from them the new-physics 
parameters sy. Most of these observables have been discussed in recent phenomenological 
analyses, usually in the framework of the type II 2HDM ]38| , |39| , [4C| , 41 1, but also in the 
type III 2HDM @. 

For that purpose, a good control of the hadronic decay parameters is necessary. These 
usually involve large theoretical uncertainties whose treatment is not well defined. In 
our work we use the statistical approach RFit [EjJ, which has been implemented in the 
CKMfitter package Q]. The new-physics parameter space is explored, assigning to each 
point the maximal relative likelihood under variation of the theoretical parameters which 
are not shown. Theoretical uncertainties are treated by defining allowed ranges within 
which the contribution of the corresponding theoretical quantity to the A% 2 is set to zero, 
while it is set to infinity outside. This treatment implies that uncertainties of this kind 
should be chosen conservatively and added linearly. 

Another related problem is the combination of different theoretical determinations of a 
hadronic quantity, which is even less well defined. We follow the prescription given in [44]. 



However, unless commented explicitly, we only take lattice results coming from numerical 
simulations with 2 + 1 flavours. For quantities concerning the light hadrons, we consider the 



determinations recommended by the Flavour Lattice Averaging Group (FLAG) [j45j, 46]. 
The obtained values are collected in table |2[ 

For f+ n (0) the only published value with 2+1 dynamical quarks is the one from 
RBC/UKQCD [|7], ||], which however fails to fulfill the FLAG standards. On the other 



hand, there is one 2-flavour result, which fulfills the FLAG criteria [49|. Although consis- 
tent with the old Leutwyler-Roos estimate [|50|1 , based on 0(p 4 ) Chiral Perturbation Theory 
(%PT), these lattice determinations are somewhat smaller than the 0(p 6 ) analytical cal- 



culations |5l], [y, 53, 54]. We take this into account and adopt the conservative range 
/£*(()) = 0.965 ± 0.010. 



To fix the values of the relevant CKM entries we only use determinations 64, 79, 80] 
which are not sensitive to the new-physics contributions. Thus, we use the V u d value 
extracted from superallowed (0 + — > + ) nuclear /3 decays and CKM unitarity to determine 
V us = A. The values of V u b and V c b = AX 2 are determined from exclusive and inclusive 
b —7- uW[ and b — > cW[ transitions, respectively, with I = e, /i. The apex (p, ff) of the 
unitarity triangle has been determined from \V u b/V c b\, A and the ratio Am B o/Am B o (see 



section |5.2|). For the top quark mass we have adopted the usual assumption that the 



Tevatron value [67] corresponds to the pole mass, but increasing its systematic error by 



1 GeV to account for the intrinsic ambiguity in the nit definition; i.e. we have taken 
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Parameter 



Value 



Comment 



IBs/ fB d 

fD a 

foj 'fD d 

Ik/ fir 

fB 



(0.242 ± 0.003 ± 0.022) GeV 
1.232 ±0.016 ±0.033 
(0.2417 ± 0.0012 ± 0.0053) GeV 
1.171 ±0.005 ±0.02 
1.192 ±0.002 ±0.013 



Our average 
Our average 
Our average 
Our average 
Our average 



m 



58 



57] 



59, B0 



fs d \l B B o/(f Bs 
Bk 



(0.266 ± 0.007 ± 0.032) GeV |§ 



B B o) 1.258 ±0.025 ±0.043 
0.732 ± 0.006 ± 0.043 



11 



[61 



\V ud \ 
A 

\V ub \ 
A 

P 



0.97425 ± 0.00022 
0.2255 ± 0.0010 
(3.8 ±0.1 ±0.4) • 10- 
0.80 ±0.01 ±0.01 
0.15 ±0.02 ±0.05 
0.38 ±0.01 ±0.06 



b -)• ulv (excl. + incl.) Q | 
b civ (excl. + incl.) §4], |65 
Our fit 
Our fit 



m u {2 GeV) 
m d (2 GeV) 
m s {2 GeV) 
m c (2 GeV) 
m b (m b ) 
fh t (m t ) 



(0.00255 GeV 



■ 0.00105/ 

(0.00504^™) GeV 

-0.025^, 
-0.035/ 
-0.07\ 
-0.il> 1 
-0.17\ 
-0.07/ 

'165.1 ±0.6 ±2.1) GeV 



(0.105; 

(1.27; 

(4.20 ; 



GeV 
GeV 
GeV 



[66] 
[66] 

m 

[66] 



Oem 

ctK2/K12 
Oem 

A T7r2/7r£2 
Oem 



-0.0070 ± 0.0018 
0.0090 ± 0.0022 
0.0016 ±0.0014 



El Si 



[7j, [7j, 74 



/f w (0) 

9h SM 



9b,SM 



1.18 ±0.04 ±0.04 
0.46 ±0.02 
0.965 ±0.010 



_p LL_j' I I' 



51, 52, 53 



-0.42112 



0.00035 
0.00018 



0.94 ±0.02 



0.07744 



- 0.00006 
0.00008 



MM 

m 

MM 



Table 2: Input values for the hadronic parameters, obtained as described in the text, 
error denotes statistical uncertainty, the second systematic/theoretical. 



The first 



m t° 6 = (173.1 ± 0.6 ± 2.1) GeV and have converted this value into the running MS mass. 
The measurements used in our analysis are listed in table ||. 

Concerning the charged-scalar mass, we will use the LEP lower bound M H ± > 78.6 GeV 
(95% CL), which does not refer to any specific Yukawa structure [jl], [81]]. This limit assumes 
only that H + decays dominantly into Uidj and l + vi. Obviously, the bound is avoided by 
a fermiophobic (inert) A2HDM with a < 1, but all our constraints would also disap- 
pear in this case. The charged scalar could still be detected through the decay mode 

—7- VF^A, provided it is kinematically allowed. Assuming a CP-conserving scalar po- 
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Observable 



Value 



Comment 



\9 b RR \r- 

Br(T - 
Br(T - 
Br(T - 
Br(B - 
Br(D - 
Br(Z? - 
Br(D s 
Br(D s 
T(K - 
r(r -> 
logC 
Br(5 - 
T(Z -» 
Br(£ - 
Br(£ - 
Am 5 o 
Am B o 
l e if| 



► fJLV T U,j,) /Bi(t ->■ ez/ r P e ) 
-7- ri/) 

* Aw) 
■> ri/) 

— >• Tf) 

Kv)/T(t ->■ 7T^) 

+ Dtv)/BR{B -> Dtv) 
bb)/T(Z — > hadrons) 

-5- ^s7)s 7 >1.6GeV 

->• X r fiV P \ 



< 0.72 (95% CL) 
(17.36 ±0.05) x 10~ 2 
(17.85 ±0.05) x 10- 2 

0.9796 ± 0.0039 
(1.73 ±0.35) x 10~ 4 
(3.82 ± 0.33) x 10~ 4 
< 1.3 x 10~ 3 (95% CL) 
(5.58 ±0.35) x 10~ 2 
(5.80 ± 0.43) x 10~ 3 

1.334 ± 0.004 
(6.50 ±0.10) x 10~ 2 
0.194 ±0.011 
0.392 ± 0.079 
0.21629 ± 0.00066 
(3.55 ± 0.26) x 10" 4 
(10.74 ±0.16) x 10~ 2 
(0.507 ±0.005) ps" 1 
(17.77 ±0.12) ps- 1 
(2.228 ±0.011) x 10 -3 



m 
m 
m 
m 



6L 87, 88 



|g§ ||, n 

III 

in@_ 

|9| 

I! 



95] 



Table 3: Measurements used in the analysis. Masses and lifetimes are taken from the PDG jp_ 



tential, OPAL finds the 95% CL constraints M H ± > 56.5 (64i 
M A < M H ± - M w ± p|. 



GeV, for 12 (15) GeV < 



4. Tree-level decays 

4.1 Lepton decays 

The pure leptonic decays I 

W couplings and of their left-handed current structure [E6L |97J, 



I'vyvi provide accurate tests of the universality of the leptonic 

. The exchange of a 



charged scalar induces an additional amplitude mediating the decay of a right-handed 
initial lepton into a right-handed final charged lepton; in standard notation 



9g, this 

scalar contribution gets parametrized through the effective low-energy coupling = 



M 



H± 



Its phenomenological effects can be isolated through the Michel parameters 
governing the decay distribution, 

3 . 1 



o, n 



2N 



Re(<? 



rr) ' 



2N 



\9rr\\ & 



3 
4 



8iV 



\9rr\ 



and in the total decay width 



T(l ->■ l' v v vi) 



G F 
192vr 3 



mf N 



V mi 



m v ( m 2 v \ 

if l 9[ — ) 

mi \mf J 



rRC 



(4.1) 



(4.2) 
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where f(x) = 1 — 8x + 8x 3 — x 4 — 12x 2 logx, g(x) = 1 + 9x — 9x 2 — x 3 + 6x(l + x) logx, 

N = 1 + i IsLl 2 and H 



1+ V ^ 7T 

2vr 1 4 



a(m;) / 25 2 \1 I" 3 m 2 m 2 



5M 2 , 



(4.3) 



Since the scalar couplings are proportional to lepton masses, the decay r — > fiv^r is the 
most sensitive one to the scalar-exchange contribution. The present bound |<?|>-rI t - s> m < 0.72 
(95% CL) H translates into \q\/M h ± < 1.96 GeV" 1 (95% CL). A better limit can be 
obtained from the ratio of the total r decay widths into the muon and electron modes. 
The universality test \g^/g e \ 2 = |Br(r -4 /i)/Br(r — > e)\\f(m 2 /m 2 )/f(m 2 ^/m 2 )\ = 1.0036± 
0.0029 ||, H implies: 



M H ± 



< 0.40 GeV" 1 (95% CL). (4.4) 



4.2 Leptonic decays of pseudoscalar mesons 

Information about new-physics parameters can be also extracted from leptonic decays of 
pseudoscalar mesons, P + — > which are very sensitive to H + exchange due to the 

helicity suppression of the SM amplitude. The total decay width is given by 2 

/ \ 2 

T{P± -4 Z+n) = GWflW^ 2 -JkU-^-j (1 + C 2 ) |1 - A,/ , (4.5) 

where i,j represent the valence quarks of the meson under consideration. The correction 

/ m p ± \ 2 ?,.m,, . + omj . 
13 \M H ±) ^ m Ui +m d] 1 ' 

encodes the new-physics information and 5^ 2 denotes the electromagnetic radiative con- 
tributions. These corrections are relevant because the additional photon lifts the helicity 
suppression of the two-body decay, thereby compensating in part for the additional elec- 
tromagnetic coupling, and the two processes are not distinguishable experimentally for 
low photon energies. Their relative importance therefore increases for decreasing lepton 
masses. 

The correction Aij is predicted to be positive in model I, negative in model X and 
can have either sign in the models II and Y, depending on the decaying meson, while it 
is of course absent in the inert scenario. In the more general A2HDM it is a complex 
number with a real part of either sign. To determine its size one needs to know | | and 
a theoretical determination of the meson decay constant. 

The SM as well as the 2HDM contribution to this class of decays start at tree level. 
Therefore they can be assumed to remain the dominant contributions, relatively indepen- 
dent of a possible high-energy completion of the theory. Electroweak loop corrections are 
of course expected and they could be sizeable in some cases, for example in supersymmetry 



at large values of tan (3 [101, 102] 



2 The normalization of the meson decay constant corresponds to f n = \/2-FV = 131 MeV. 
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4.2.1 Heavy pseudoscalar mesons 



E 

-0.05 




o.6 grt- 



0.1 




0.2 0.4 0.6 



Refi^/M*) 



Figure 1: Constraints in the complex ?*Su,<i/Afjy± planes from B — > rf (1e/tj and D [iv (right), 
in units of GeV~ 2 . The colour code indicates confidence levels (1 — CL). 



The leptonic decays of heavy pseudoscalar mesons that have been measured up to now 
are B — > tv, D s — > jj,v, D s — > tv and D — > [xv. The radiative corrections for the leptonic 



decays of heavy mesons have been estimated in [103], and are already taken into account 
in the experimental values given in table |3|; therefore the electromagnetic correction is set 
to zero in Eq. (4.5). 

In B and D decays the function Ay can be approximated by neglecting the contribution 
proportional to the light quark mass, because m u /nib < md/rn c ~ 0(1O~ 3 ). Therefore the 
relations 



A 



B 



A 



m 



cd 



D 



(4.7) 



hold, leading to a direct constraint on these combinations. While for — > tv the helicity 
suppression is absent, the corresponding phase space is small and there are two neutrinos 
in the final state, which is why D — > tv has not been measured up to now. Nevertheless, 
the upper limit set by CLEO |84| starts to become relevant in constraining our parameters: 
|1 — A c d\ < 1.19 (95% CL). The present experimental limit on B — > /j,v gives |1— A u },\ < 2.04 
(95% CL). The information obtained from the decays B — > tv and D — > fiv is shown 
in figure [j]. The broad dark red (black) ring in the middle reflects the fact, that the 
systematic error is dominant in these constraints, leading to a large amount of degeneracy 
for the 'best fit value'. To infer a limit at a certain confidence level, the corresponding 
number of rings has to be included, for example for 95% up to the yellow (light grey) 
corresponding to 1 — CL = 0.05. The resulting 95% CL constraints, |1 — A u ^| £ [0.8,2.0] 
and |1 — A C( i\ € [0.87, 1.12], translate into allowed circular bands in the tfs U: d/Mjj ± complex 
planes. For real Yukawa couplings there is a two-fold sign ambiguity generating two possible 
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Re(C u Ci/Mg) 

Figure 2: 95% CL constraints in the complex <;*<; u /Mjt± plane from D s — > (r,fj,)v, in units of 
GeV~ 2 , using B — > rv to constrain <^*q;/M^ ± . 




0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 



Re(#&)/M| Re(ff&)/M| 

Figure 3: Constraints from D s — > tv t (left) and D s —> fiv^ (right), in units of GeV~ 2 , under the 
assumption of real parameters The grey bands correspond to 95% CL. Also shown are the cuts 
for the 2HDM of type I/X (dashed line) and II (lighter grey area, tan (3 G [0.1, 60]^. Finally, the 
four black regions are the possible allowed areas considering the information coming from B — > rv T . 



solutions, the expected one around Aij = (the SM amplitude dominates) and its mirror 
around Ajj = 2, corresponding to a new-physics contribution twice as large as the SM 
one and of opposite sign. The real solutions are q*q d /M^ ± G [-0.036,0.008] GeV -2 and 
[0.064,0.108] GeV~ 2 , and ^q u /M 2 H± G [-0.037,0.037] GeV~ 2 and [0.535,0.609] GeV" 2 . 

In D s decays we get |1 - A cs | G [0.97, 1.18] from D s ->• fiv and |1 - A cs | G [0.98, 1.16] 
from D s —> tv. Here the situation is a bit more complex, because m s /m c ~ 10% and 
the light-quark term in the A cs function cannot be neglected since this suppression could 
be compensated by the different <jj. Therefore there is no direct constraint, neither on 
<^i<; u /Mjj ± nor on ^*<^/M^ ± , only a correlation among them. For that reason, we use the 
additional information from B — > tv to constrain the parameters which are not shown. 



This suffices to render the influence of the mass-suppressed term subdominant. 

If CP symmetry were only broken by the CKM phase, the parameters cry would be 
real. In this case, the constraints from D s — > tv t and D s — > fj,Vu can be visualized as 
shown in figure |3], plotting the correlation between the two real parameters. The two grey 
bands are associated with the two possible solutions around A cs = and A cs = 2. The 
different models with ^-symmetry correspond to cuts in these plots. The plots show the 
small influence of the term proportional to the strange quark mass, as long as the couplings 
are of the same order. Using the constraints on ^*^/M^- ± from B — > tu, one finds for 
the other coupling combination the two real solutions ^*? U /M^ ± G [—0.005,0.041] GeV -2 
and [0.511,0.557] GeV -2 , at 95% CL, which agree with the corresponding constraints from 
D — > [iv. Putting together all the information from leptonic B, D and D s decays, the real 
solutions are: 

[-0.036, 0.008] GeV" 2 , 
[0.064,0.108] GeV -2 , 



[-0.006,0.037] GeV" 2 
[0.511,0.535] GeV -2 . 



4.2.2 Light pseudoscalar mesons 



(4.8) 



r 1 (jj 




kJ! o 
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era : 
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o 1. 



r(K^nv)/ 
r(;i— >(iv) 



r(x->Kv)/ 
r(t— >7tv) 



Re(CCi/Mfi) 



Re(Q{ d )/M 2 H 

Figure 4: Constraints in the complex plane (s* Sd) / Mfj± , in units of GeV~ 2 . Left: Full regions 
allowed at 95% CL for K/ir — > \xv (upper plot) and r — » K/ttv (lower plot). Right: 95% CL 
constraints in the interesting region (from the global fit) for both constraints, using D — » \w to 
constrain s*<r„/Mj^ ± . 



Due to the cancellation of common uncertainties, lattice calculations of the ratio fx/ fn 
are more precise than the determinations of the individual decay constants. This ratio can 
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be extracted experimentally from two different ratios of decay widths: 



T{K -> nu) 
T(t -> Ku) 



m K 



m T . 



1 



m 2 /m 2 K 



1 



m 



- ml/ml 

I 2 \ 2 



V„ 



V u d 



1 - A, 



A' 



r,, 



(i + &^ 2 ) 



l 



i - 



A u d 

2 



1 - A 



AtK2/K12) rKl2/irt2 
Oem i Oem 



ud 



j7T7r2/7r£2 



(4.9) 



(4.10) 



0.0004± 



r(r ->■ 7r^) ^l-m^/m^y F ud 

where o em is given m table y and o cm 
0.0054. 

The new-physics corrections are dominated by A us ~ Qin^ / 'M^± . As m 2 K /m 2 B ~ 
1%, the scalar contributions to these decays are much smaller than for the heavy mesons. 
However, the good experimental precision achieved provides interesting constraints, as 
shown in figure ||, which are dominated by the K^/ir^ ratio. At 95% CL, one finds 
|1 - A us | G [0.984, 1.017] from and |1 - A us \ G [0.965, 1.025] from the ratio r -> 

uK/tt. The real solutions are then, ?*^/M|- ± G [-0.07,0.07] GeV~ 2 or [8.14,8.28] GeV~ 2 . 
The larger real solution is already excluded by the B — )• rv data. 

4.3 Semileptonic decays of pseudoscalar mesons 

Semileptonic decays receive contributions from a charged scalar as well, but in this case the 
leading SM amplitude is not helicity suppressed, therefore the relative influence is smaller. 
In addition, there are momentum-dependent form factors involved. The decay amplitude 
M —> M'Wi is characterized by two form factors, f+(t) and fo(t) associated with the 
P-wave and S-wave projections of the crossed-channel matrix element (0\ui , y^dj\MM'). 
The scalar-exchange amplitude only contributes to the scalar form factor; it amounts to a 
multiplicative correction 

fo(t) = f (t) (l + <%i) i ( 4 - n ) 



where 



in,, 



(4.12) 



The determination of the CKM matrix element \Vij\ is not contaminated by the new- 
physics contribution, because it is governed by the vector form factor. One measures 
the electron mode M — > M'ev e , where the scalar contribution is heavily suppressed by the 
electron mass, determining the product |V^| |/+(io)|, with to = for light-quark transitions 
and to = (mM — niM') 2 for heavy quarks. A theoretical calculation of |/+(io)| is then 
needed to extract \Vij\. The sensitivity to the scalar contribution can only be achieved 
in semileptonic decays into heavier leptons. Whenever available, one can make use of 
the differential decay distribution to separate the scalar and vector amplitudes. In any 
case, theoretical determinations of the scalar and vector form factors are needed to extract 
information on Sij. 

4.3.1 B -> Dtv t 

To reduce the uncertainty from the vector form factor, let us consider the ratio 



Br(B -»• Dtv t ) 
Bt{B ->• Dev e ) 



a Q + ai (m| - m 2 D ) Re(<5 cb ) + a 2 (m| - m 2 D ) \8, 



cb 



(4.13) 
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The coefficients Gtj, which contain the dependence on the strong-interaction dynamics, have 
been studied recently and parametrized in terms of the vector form-factor slope p 2 and the 
scalar density A(vb • vd) = A, assumed to be constant [39, 104]. We make use of these 
parametrizations, taking for the two parameters the values indicated in table ||. The 
function A(vb • vd) ot /o(t)//+(t) has been studied in the lattice, in the range vb-vd = 1- 
1.2, and found to be consistent with a constant value A = 0.46 ± 0.02, very close to its 
static-limit approximation (tub — mo) /(tub + mo) J75[. 




ReK d ?;/Mg)/GeV-' 
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Figure 5: Constraints from B — > Dtv t , in units of GeV~ 2 , plotted in the complex plane for 
S*Sd/Mff± (left) and <;*<; u /Mjj± (right), using D — > \w and B — > tv to constrain the combination 
not shown, respectively. The colours indicate 1 — CL, the red lines the constraint (95% CL) for 



We obtain once more a correlation between s*s u /Mfj± and ^*^d/M%±, where the term 
proportional to the charm quark mass is in general potentially more important than in the 
type II model. The results are shown in figure |5| for both parameter combinations. As can 
be seen there, the constraint on <^*/Mjl ± is consistent with the information coming from 
B — > tv and the leptonic decays of light mesons, but does not constrain this combination 
further as long as only the information of B — > Dtv t is used. The red lines indicating the 
constraint for g*^ — > 0, however, show that the semileptonic decay can exclude a small 
region around (0.08,0), once that combination is bound to be small. We will use this to 
exclude the second real solution for ?*<^/Mj^ ± with aid of the processes e^,Z — > bb and 
r — > [ivv (see figure Also, when plotted in the complex ^/ML plane, it becomes 
apparent that this constraint is important to exclude the second real solution allowed by 
Di s \ — > tv decays, already using only the information from leptonic decays in addition (see 
again figure ||) . 

Considering the limit of real ?/'s, the correlation between the real parts is visualized in 
figure [f| together with the cuts corresponding to the different models with Z2 symmetries. 
The plot shows that the mj, and m c terms have potentially similar influence in this case. 

It has been pointed out in [105] that measuring the spectrum instead of just the 
branching ratio will increase the sensitivity of this channel. This, however, has not been 
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done up to now, due to lack of statistics. 




Re(£*£/)/Ml 

Figure 6: Allowed regions for Re{<;* / MZj± and i?e(^*^ u )/M| r± from the process B — > Dtv at 
95% CL (grey), in units of GeV~ 2 , assuming that their imaginary parts are zero. The projections 
for the 2HDMs of types I/X ( dashed line ) and II (lighter grey area, tan /? 6 [0.1, 60],) are also shown. 



4.3.2 K iriu 

In semileptonic kaon decays the Callan-Treiman theorem 1 106 , 107] allows to relate the 
scalar form factor at the kinematic point t CT = m 2 K — m 2 to the ratio of kaon and pion 



decay constants: C = /o(t CT )//+(0) = f+(o) ^cn where A 
small xPT correction of 0[m 2 / (An f^) 2 



CT 
1 



CT 



-3.5 ± 



10 



3 is a 



108 , |54] , 109 1 . Using a twice-subtracted dispersion 



relation for fo(t) |110| , the constant C has been determined from the data by KLOE 
|91[, KTeV 



and NA48 [111]. In the average quoted in table || the NA48 result has been 
excluded because it disagrees with the other two measurements by more than 2a. 

In the presence of charged-scalar contributions, the scalar form factor gets modified 



as indicated in Eq. ( 4.11 ), inducing a corresponding change in C. Taking into account 
that the analyzed experimental distribution is only sensitive to |/o(t)| 2 , to first order in 
the new-physics correction 5 US , the measured value of C corresponds to 



bgc = log (^fj^) + AcT ) + Re 



(4.14) 



The resulting constraint on the real part of c d ^*/M^ ± is shown in figure [7|, leading to 



Re 



-0.16,0.30] GeV" 2 (95% CL) 



(4.15) 



which is in agreement with the previous constraints, but with larger uncertainties. This 
might change in the near future, due to improved lattice determinations of /+(0) and 
fx/f-K, as well as improved experimental precision, e.g. from NA62 or KLOE-2. 

4.4 Global fit to leptonic and semileptonic decays 

Combining the information from all leptonic and semileptonic decays discussed before, 
one gets the constraints shown in figure ||. |<^*/M^ ± | 



is bounded to be smaller than 
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Figure 7: Constraint from the direct measurement of log C , in units of GeV 



~ 0.1 GeV -2 (95% CL) from these decays alone, while for q u q^ /M^ ± the constraints are 
relatively weak, due to the similar masses of the mesons in the leptonic decays. Note 
that in both cases there are two real solutions. For the combination ? U <^*/M|^ ± , one real 
solution is excluded in the global fit at 95% CL, while the other, including the SM point 
of vanishing couplings remains allowed. As mentioned before, this exclusion is due to 
B — > Dtu in combination with the constraint on ^^*/M^- ± . For the latter, the situation 
is more complicated. The second solution remains allowed, due to the overlapping of the 
two main constraints in both regions and the weak constraint on /M^ ± derived from 
semileptonic decays. However, using in addition the information coming from leptonic r 
decays in (4.4), the lower Higgs mass bound from LEP and the constraint from €k, Z — > bb 



(see section |5.1| ) in a conservative way, \<i u q*\/M'r I ± < 0.01 GeV 2 , the second real solution 
for SdS* /Mjj ± is excluded as well by B — > Dtu. 




! D— >\iv 




R(B^DIv)(+B^tv) 



-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Re^/MS 



Figure 8: /M"^ ± (left) and ?„<j* / M^± (right) in the complex plane, in units of GeV~ 2 , con- 
strained by leptonic and semileptonic decays. The inner yellow area shows the allowed region at 
95% CL, in the case of ^* /M^± using additional information (see text). 
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5. Loop-induced processes 



For processes where new-physics contributions appear only through quantum loop effects, 
the situation becomes obviously more difficult, regarding not only the calculation but also 
the interpretation of the results. If the SM amplitude is also mediated only by loops, the 
relative importance of the charged-scalar contributions is expected to be higher, but this 
implies also a higher sensitivity to the framework in which the A2HDM is eventually to 
be embedded in. In the following we make the assumption that for the observables under 
discussion the dominant new-physics corrections are those generated by the charged scalar. 
Moreover, since no significant signal for new physics has been found up to now in flavour 
observables, we assume these effects to be subleading with respect to the SM contribution. 



5.1 Z -> bb 



The high-precision data collected at LEP and SLD has made it possible to accurately test 
the SM electroweak loop corrections at the Z scale, providing information on the Higgs mass 
and useful constraints on many new-physics scenarios. While most Z-peak observables 
are only sensitive to the gauge-boson selfenergies, the decay Z — > bb provides valuable 
information on fermionic vertex corrections induced by charged-current exchanges. Since 
Vfb ~ 1, those loop diagrams involving virtual top quarks generate quantum corrections 
to the Zbb vertex, which are absent in the Zdd and Zss vertices. These corrections are 
enhanced by a factor mf, allowing for a quite accurate determination of the top quark 
mass |112| , 113]. The same arguments apply to the charged-scalar contributions present 



in the A2HDM, providing a sensitive probe of the corresponding H + tb coupling. For very 
large values of |?^| this decay would also be sensitive to contributions from neutral scalars 



114]; we don't consider this possibility here. However, given a not too small value for q, 



(semi-)leptonic decays can be used to exclude that possibility. 

Therefore, we assume the dominance of charged-scalar effects in the following, allowing 
only for |c^| < 50. We disregard the information coming from the forward-backward polar- 
ization asymmetry A^, because the scalar-exchange contributions to A^ are small compared 
to the present uncertainties. 

It is convenient to normalize the Z —> bb decay width to the total hadronic width of 
the Z, because many QCD and electroweak corrections cancel in the ratio, amplifying the 



sensitivity to the wanted vertex contribution [112]. Within the A2HDM, this ratio can be 



written as |TJ, HI' llZl 



Rb = 



T(Z 



T(Z — > hadrons) 



1 + * q QCD 

Sb 



(5.1) 



where 



», = [(St - S?f + (it + sf ) 21 



Si, 



g^b,t 



q > 



(5.2) 
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Figure 9: Constraint from Rb in the \<> U \ — M H ± plane (M H ± in GeV units), allowing for |<^| < 50. 



with C 6 QCD = 1.0086 being a factor including QCD corrections |116| . The A2HDM contri- 
butions are encoded through the effective left- and right-handed Z66-couplings: 



9b 



9b,SM + 



9b, SM 



w 



16vr 2 
V2G F M^ 



16vr 2 



m t 

m b 
M 2 W 



fi(t h ) + ^f 2 (th) 

fl(t h ) + ^f2(t h ) 



(5.3) 
(5.4) 



where t h = mf/M^, fi(t h ) = [t 2 h -t h -t h \ogt h ]/(l-t h ) 2 and the function f 2 (t h ) governing 
the NLO correction is given in j77|. If running quark masses fht(Mz) and fhj,(Mz) are used, 
this NLO QCD correction is small. The light-quark coupling contribution Sf, = 1.3214 



-L,R 



1 76, 77] is not sensitive to the new-physics effects. The SM values of the couplings g b ' SM , 



given in table have been computed removing the Z — > bb information from the standard 
electroweak fit [76, 77]. 

In contrast to the leptonic and semileptonic constraints discussed before, here the 
parameters enter directly, allowing to bound them without information on \q\. The 

constraint resulting from the input values in tables |2| and ||| is shown in figure ||. The 
constraint is plotted in the |? u | — M H ± plane, as obviously it is much weaker for |<^|, due 
to the relative factor irib/mt. For large scalar masses, the constraint weakens as the effects 
start to decouple, reflected in limt ;i ->o fxpith) = 0- I n the range of scalar masses considered, 
it leads to a 95% CL upper bound \q u \ < 0.91 (1.91), for M H ± = 80 (500) GeV. The upper 
bound increases linearly with M H ±, implying 



< 0.0024 GeV" 1 + 



0.72 



< 0.011 GeV 



-1 



(5.5) 



M H ± M H ± 

where we have used the lower bound on the charged-scalar mass from LEP searches, M H ± > 
78.6 GeV (95% CL) |, 0. Combined with the limit on \q/M H ± \ from leptonic r decays, 
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this already constrains the combination |^ n <^*|/M^- ± much stronger than the global fit to 
(semi)leptonic decays, leading to 



M 2 



< 0.005 GeV -2 , (5.6) 



however only with the additional assumptions of |c^| < 50 and charged-scalar effects dom- 
inating the new-physics contributions to Rb- The range allowed for |<^| in the fit does not 
influence the upper bound on \q u \, apart from the exclusion of neutral-scalar effects, since 
both contributions can only lower the value for Rb and both are allowed to vanish in the 
fit. Therefore the upper limit stems from points with |c^| = 0. 

5.2 B°-B° mixing 

The mixing of neutral B mesons is very sensitive to charged-scalar effects, as the leading 
contribution stems from top-quark loops, rendering the new-physics and SM contributions 
comparable. Besides the high precision of the measurement for the mass difference Am^o, 
the B® mixing is especially interesting due to the observed tension in its phase |65]| . 
In the usual 2HDMs with a Z 2 symmetry the scalar couplings are necessarily real, leading 
to a vanishing contribution to this phase. However, the complex Yukawa couplings ? Uj d 
of the A2HDM provide a potential new-physics contribution, which could account for the 
experimentally observed phase. 

In the SM, the calculation is simplified by the fact that only one operator contributes, 
denoted VLL below. In the presence of a charged scalar, an enlarged effective Hamiltonian 



has to be considered, involving a basis of eight operators [117, 118, |119| , 12C 1 : 
VLL,VRR = ^V L , R b«) {dP-fVL^fi) 

o\ K = {d a ltl V L b a ) (&-fP R tf*) , 
Q\ K = (d a V L b a ) (dPv R bP\ , (5.8) 

SLL,SRR = (d a V L , R b a ) (d^V L ,Rl 



sll,srr = ^ a a^V L , R b a ) (dP<T^V L ,R 



with a, /3 being colour indices and a^ v = ofT^T^]- We have written the effective Hamil- 
tonian relevant for B®-B® mixing; the mixing of B® mesons is described by the analogous 
expression, changing the label d to s everywhere. 

We have performed the matching of the underlying A2HDM and the low-energy effec- 
tive Hamiltonian at the scale \xtw ~ M\y,mt. The resulting Wilson coefficients, given in 
the appendix, reproduce the SM result as well as the matching for the 2HDM in the limit 



ma — > 0, given in [121]. As noted above, the contribution of the A2HDM to Cyll(iMw) 
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is an 0(1) effect. For that reason, we calculate this contribution at NLO, implementing 
the results of (l2| within the A2HDM. 3 O wing to their chirality structure, the remaining 
Wilson coefficients are all suppressed by powers of the light-quark mass (m s in the 
B® case), except Cg RR which is proportional to m%. Restricting the parameter ranges to 
\s u \ £ [0,5] and |^| G [0,50], the ratio \Ci{ntw) I 'Cvll(^w)! is then below two percent for 
all operators apart from 0^ RR . Since the matrix elements for the B° mixing do not contain 
the large (chiral) enhancement factors present in the kaon system, this allows us to restrict 
ourselves to two operators only. Moreover, the ratio Cg RR /CvLL is a small quantity (10% 
at most for < 25, still below 40% for |q| = 50) and therefore a leading-order estimate of 
the contribution is enough for our purposes, while the dominant VLL contribution 

is included at NLO. 

The strong (m s — m^jMy/ suppression of SU(3)-breaking effects implies that, for 
the parameter ranges considered, the ratio Am B o/ Am B o is unaffected by charged-scalar 

s d 

contributions and can be used in the CKM fit. Note, however, that in the limit |^| ^> 
50, <C 1, which corresponds to the large-tan/? scenario in the type II model, the 
contribution from VRR might become the dominant new-physics correction to B® mixing, 
but remains small compared to the SM one. 

We use the ratio Am B o/Am 8 o to determine the apex (p, ff) of the unitarity triangle, 
and bound the charged-scalar parameters with the B® mixing information. The resulting 
constraint from Amgo in the Mjj±— |? n | plane is shown in figure 10, using the scales 
fj-tw = m t and Hb = 4.2 GeV. The error includes the variations in the CKM parameters, 
f B o, B B o and the experimental uncertainty. The leading VLL contribution depends on 
|? u | 2 only, while Cg RR is proportional to = \s u \ |^d|e l¥> , (p being the relative phase 



between the two Yukawa couplings. To determine the allowed region shown in figure 10 
we have varied ?d in the range |q| < 50 and <p> G [0, 2tt]. 

Interestingly, the dominant contribution to a possible phase shift in the mixing is 
also the one from The factor My^Do(mt, M H ±) (see appendix) varies between zero 

and ~ —3% for scalar masses between 50 and 500 GeV, while Am^mf/M^ ~ 10%. For 
relatively large values of the product |?„<f d | (> 20) this factor can contribute sizeably to 
the B° mixing phase, as long as M H ± is relatively small. The sign of the shift is obviously 
not fixed, but depends on the sign of the relative Yukawa phase ip. As long as \^d\ is not 
too large, the effect is the same in B^ and B®. 

The DO experiment has measured very recently |3(| a like-sign dimuon charge asym- 



metry leading to A b , = -0.00957 ± 0.00251 ± 0.00146, which differs by over three standard 



deviations from the SM prediction |122| , 123]. The measurement includes contributions 



from B° d and B° s mesons, corresponding to A b sl = (0.506 ± 0.043) a d sl + (0.494 ± 0.043) a s s 



with (q = d, s) 



. r g \ \r q i Ar B o 

a q , = Im -Jf = J-^- sm(b a = tand> a , (5.9) 

\M?J \M? 2 \ Vq AM B o V<? ' 1 > 



where Mf 2 — ^T q l2 = (Bq\"H^ 2 \B®). While this result needs certainly confirmation, we 
will explore some of its consequences for the parameters of the A2HDM in the follow- 

3 Note, that there are several smaller errors in that paper, most of which have been pointed out in [p8|. 
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Figure 10: The 95% CL constraint coming from Am B o in the M H ± - \<; u \ plane for |q| G [0,50], 
varying in addition the relative phase if in [0, 271"]. The excluded area lies above the dark (red) region 
only. In yellow the allowed area for ^ = is shown. 



ing. Using the current experimental value for the asymmetry in the B® system, a*?, = 
-0.0047 ± 0.0046 §5j, the measured value of AM B o and the SM prediction for AT B o, the 
DO asymmetry implies sin^ s = —2.7 ± 1.4 ± 1.6, showing that the central value of this 
measurement is incompatible with the assumption of negligible influence of new physics 
on rj 2 , while the uncertainties are large enough to allow every value for the mixing phase 
at 2a. Using in addition the direct measurement of a s sl through B® — > fi + D~X decays by 
DO pi , = -0-0017 ± 0.0091, results in sinc/> s = -1.7 ± 1.1 ± 1.0. Note that (part 
of) the observed deviation may also be due to the possibility of bad convergence of the 



operator product expansion (OPE) |125|, 126], related to the relatively low effective energy 



scale mj) — 2m c . However, no signs for a breakdown were found in the above calculation. 
Note also that in [127] it has been argued that such a large value violates a "coherence 
bound" derived by demanding monotonicity of the Stokes vector in the B s system. The 
possibility of NP influence on the rate as an explanation for this measurement has recently 



been discussed in |128| , |129| , |130f| . The authors of |128| conclude, that most of the possible 



operators are strongly constrained by other processes (including the one discussed in |129f| ), 
leaving little space for an 0(1) contribution to rf 2 - 

Hints of a large (j> 8 value have been also obtained previously from B® — > J/ip(j) decays 
1 131, 132, 133], where the extraction of the phase might however be influenced by contri- 
butions to the decay amplitude: in the SM, one of the reasons why this decay is "golden" 
is the fact, that the potentially relatively large penguin contributions have the same phase 
as the leading (colour-suppressed) tree amplitude, and therefore do not spoil the extrac- 
tion of the mixing phase from the time-dependent CP asymmetry. However, this is no 
longer true in the A2HDM: the charged-scalar penguin contributions include terms similar 
to their leading SM counterparts, with an additional factor of ?*^m;,mt/M^ ± ~ 0(1), 
thereby providing a second weak phase in the decay amplitude. Quantitatively assessing 
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the influence of these contributions would require a reliable calculation of the correspond- 
ing matrix elements, which is however not available; we are thus left with the possibility 



of a semi-quantitive analysis only, e.g. along the lines of [134], which we however do not 
consider here. 



The SM predicts a very small positive value for . 
for (f>d- The theoretical values quoted in [122] are 4- 



and a much larger and negative result 
= 0.24° ± 0.08° and <j) d = -5.2° ±lf a . 



Sto(«/(|A I |StoC$ M )) 
150 



Sin((Sj)/(|A rf |Sin(^")) 





ArgffiM/rad 



Figure 11: Dependence of sin (j> q /(\A q \ sin<^ M ) on ip = arg (<^ d ). 



Assuming that the charged-scalar contributions are the only relevant new-physics ef- 
fects, we can analyze the possibility to accommodate a large (p s phase within the A2HDM. 
In figure 11 we plot the allowed range for sin A 9 | sin<^ M ), where A q = Mf 2 /Mf^ M , as 
a function of the relative Yukawa phase <p = arg (<^? d ). The other scalar parameters have 
been varied in the ranges |q| G [0,50], M H ± G [80,500] GeV, and according to the 
allowed range from £k,Z — ► bb, which includes only values for (\s u \,M H ±) which lead to 
acceptable values for Am Sj( i- While it is indeed possible to obtain a large value of cj> s , the 
predicted equality of A s and implies a strong anti-correlation of sin^/(|Ad| sm</>^ ) 
and sin <j> a /{\ A s | sin(^ M ), due to the different sign (and size) of </>^ M and <ft^ M - This leads to 
a prediction for the sign of ai, which could be verified/falsified, once higher experimental 
precision is achieved. As can be seen, the preferred negative sign for the a s sl asymmetry 
implies ip £ [7r/2,7r], [Sir/2, 2-ir], and for possible large values the Yukawa phase should not 
be close to 0,7r (obviously). 

Figure |l2| shows the dependence of sin S /(| A s | sin<?!>g M ) with |^| (left) and M H ± 
(right), varying the remaining parameters within their allowed ranges. If large values for 
the a s sl asymmetry are confirmed (within the physical range | sin<^ s | < 1), this would point 
towards large values of |q| and small charged scalar masses. Finally we show in figure 13 
the plots from figure 11 again, restricting the product < 20 (see section 5^). The 

corresponding maximal asymmetry is correspondingly smaller, but still relative factors up 
to ~ 60 are allowed for B s with respect to the SM. 

Additional contributions to (j) s could be induced by neutral scalar exchanges, through 
the effective FCNC operator in Eq. fl2.13j ) appearing at the one-loop level. Also, a sizable 
Yukawa phase ip = arg(?*? d ) could generate observable signals in other CP- violating ob- 
servables not yet included in our analysis. A detailed discussion of these effects and their 
corresponding constraints on the model parameters is postponed to future work. 
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5.3 K°-K° mixing: e K 

The AS = 2 effective Hamiltonian is described by the same basis of four-quark operators 
given in ( |5.8| ), changing the flavour b to s everywhere. However, the small light-quark 
masses rrid and m s suppress now the contributions from all operators except C VLL . Another 
difference with respect to B° mixing is that, owing to the different CKM factors, one needs 
to consider the virtual contributions from top and charm quark exchanges within the box 
diagrams: 

C 2 M 2 

n? s S=2 = -^f {\ 2 C% L M + \lC$lM + 2\ t \ c C$ L M} (d 7 ,V L s) {SrfVLs) ■ 

(5 " 10) 

Since At = V t * d V ts ~ A 2 X 5 while Xt = V* d V cs ~ A, in spite of the m 2 /m 2 relative suppression, 
the charm loop gives the dominant short-distance contribution to Anix- There are in 
addition large corrections from long-distance physics, which make it difficult to extract 
from Attik useful constraints on the new-physics amplitude. 

More interesting is the CP- violating parameter ex, which can be written in the form 

K e e**« Im(Mi 2 ) 

e K = 7= 7 , (5.11) 

v/2 Am K y 1 



where K e = 0.94 ± 0.02 takes into account small long-distance corrections |135| , |7§| 1. The 
top and charm contributions are now weighted by less hierarchical CKM factors Im(Af) ~ 
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A 4 Im(A c Af) ~ A 4 Im(A 2 ); the mass hierarchy compensates for this, implying that the top 
quark gives the most important contribution to ex- 

The relevant Wilson coefficients Cy LL , containing the SM and new-physics contri- 
butions, are given in the appendix. The corrections induced by the charged scalar are 
proportional to |? M | 2 and |?«| 4 . All contributions from the coupling q are absent in the 
limit rrid, s = 0. The matrix element (K°\H^ = ' 2 \K ) is parametrized through the hadronic 
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Figure 14: 95% CL constraints from ck- 

quantity f^B^-- We use the numerical values of fx/fw and Bk, given in table ||, together 

with the phenomenological determination of the pion decay constant from r(-7r + — > 

/tt = 130.4 ± 0.04 ± 0.2 MeV [S6|. Figure 14 shows the constraint obtained from €k in the 

plane M H ± — \q u \. It is very similar to the one extracted from Z — > bb, and even slightly 

stronger. 



5.4 B X s 7 



The radiative decay B — > X s j has been calculated at NNLO in the SM, leading to the 
prediction Br(B -> X s -/) S m = (3.15±0.23) x 10" 4 p6| . In the 2HDM the decay amplitude 
is known at NLO 1137] , |138| , |77| , |139| . Following the steps given in |140(| , one can express 
the branching ratio as 



Br(5 -> X sl ) Ey>Eo = Br(B -> X c eu) 



cxp 



VtsVtb 



Vcb 



6a 



[P(E ) + N(Eq)] , (5.12) 



where the phase-space factor C B = \V ub /V cb \ 2 T{B -> X c eu)/Y{B -> X u ev) = 0.580±0.016 



1 141] accounts for the m c dependence of Br(£? — > X c eu). Normalizing the result with the 
B — > X c ev transition, cancels the leading non-perturbative corrections of order A 2 /m 2 
and minimizes many sources of uncertainties, such as those generated by the CKM quark- 
mixing factors, the dependence on m\ and the sensitivity to m c . The subleading non- 
perturbative contributions are contained in N(Eq), which includes corrections of 0(k 2 /m 2 ) 
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IIH], C>(A 3 /m|), 0(A 3 /m b m 2 c ) g|| and 0(a s A 2 /(m b - 2E ) 2 ) pj . The relevant com- 



bination of CKM factors is given by 



2 



1 + X 2 (2p- 1) + X 4 (p 2 +fj 2 - A 2 ) + 0(X 6 ) = 0.963 ±0.002 ±0.005, (5.13) 



where the sensitivity to the apex (p, rj) of the unitarity triangle is suppressed by two powers 
of A. 

For m s = the effective low-energy operator basis remains the same as in the SM. The 
modifications induced by new-physics contributions appear only in the Wilson coefficients, 
which are included in the perturbative part P(Eq): 

Cf i^w) = Ci,SM + k«P Ci,uu - (<£$d ) Ci,ud , 

where = |? u ||Qi|e l</3 , V 9 being the relative phase. The virtual top-quark contributions 

dominate the coefficients Ci iUU and C% u d\ their explicit expressions as a function of mt can 
be found in [77]. Depending on the value of the phase ip, the combined effect of the two 



terms Ci )UU and Ci tU d can be rather different. For instance, these two terms tend to cancel 
each other in the type I model where ip = 0, while in the type II version with (p = ir they 
add constructively. 

Since the new-physics contribution is only calculated up to NLO, terms in the branch- 
ing ratio of 0(a 2 ) coming from the square of the 2HDM amplitude are neglected con- 
sistently. In some regions of the parameter space, leading to large new-physics effects of 
opposite sign to the SM amplitude, the cancellations between the two contributions en- 
hance the sensitivity to higher-order QCD corrections, generating in some cases unphysical 



results (for instance in the type I model at small values of tan/3) [138]. Fortunately, the 
most problematic region (large values of \s u \) is already excluded by the constraints from 
Z — > bb and Am 6 o. The inclusion of the SM NNLO contributions substantially improves 
the reliability of the theoretical predictions. 

To extract the information on the A2HDM couplings, we take into account the latest 



experimental values, given in table and use the same renormalization scales as in [140] 
(po = 160 GeV, p b = 2.5 GeV and p c = 1.5 GeV as central values and the same ranges of 
variation). We follow again the RFit approach, adding the theoretical uncertainty linearly 



to the systematic error. The resulting constraints on \<; u \ and \$d\ are shown in figure |15 
varying the charged-scalar mass in the range M H ± E [80, 500] GeV. The white areas are 
excluded at 95% CL. In the left plot, the phase (p has been scanned in the whole range 
from to 2tt; the resulting constraints are not very strong because a destructive interference 
between the two terms in ( |5,14 ) can be adjusted through the relative phase. In the range 



< 2, one finds roughly |<^||? u | < 20 (95% CL). More stringent bounds are obtained at 
fixed values of the relative phase. This is shown in the right plot, where ^ and q have 
been assumed to be real (i.e. (p = or tt). In that case, couplings of different sign are 
excluded, except at very small values, while a broad region of large equal-sign couplings is 
allowed, reflecting again the possibility of a destructive interference. 



The sensitivity to the charged-scalar mass is illustrated in figure 16, which shows the 



constraints on |<^| versus M H ± for fixed values of s u = 0.5 (left) and ? M = 1.5 (right). Again, 
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Figure 15: Constraints on ? u and qi from B — > X s j, taking M H ± € [80,500] GeV. The white 
areas are excluded at 95% CL. The black line corresponds to the upper limit from epc,Z — > bb on 
|?u|- In the left panel, the relative phase has been varied in the range ip € [0,27r]. The right panel 
assumes real couplings. 



in the upper plots the relative phase has been varied in the whole range ip E [0, 2ir], while 
the lower plots assume real couplings. Figure [17] shows the constraints on the — M H ± 



plane, for q = 0. Finally, in figure |18| we show the constraints obtained for fixed values 



of the charged-scalar mass, assuming q u and <id to be real. We reproduce in this case the 



qualitative behaviour obtained in [ 42| 



We observe that for small values of \q u \ no constraint on is obtained, because in the 
limit |<j n | —> the SM is recovered, which is compatible with the data. With growing \s u \ a 
bound on |<^| emerges, corresponding to \s u Sd\ ^5 20. For = on the other hand, a limit of 
k«l ^5 3 can be observed for large scalar masses around 500 GeV, strengthening to \q u \ < 1.3 
for smaller values of M^± . The overall constraint is relatively weak compared to the strong 
bound on M^± obtained in the type II 2HDM, due to the correlation = — 1. However, 
it can be seen from the plots with vanishing phase and/or a fixed value for that this 
strength is recovered, once some parameters are constrained independently. Comparing the 
plots with complex input parameters to their real counterparts, we observe that the effect 
of the relative phase is mainly to extend the allowed bands in a way that the excluded 
space between them is rendered allowed, too. 

We have also analyzed the CP rate asymmetry, defined as 

= BR(B -> X sl ) - BR(B -> X sl ) 
a ° P BR(B -> X sl ) + BR{B Xsj) ' 

which is predicted to be tiny in the SM. Once the constraints from the branching ratio are 
implemented in the A2HDM, the predicted asymmetry is smaller than the present exper- 
imental bounds. Thus, one does not obtain further constraints on the model parameters. 
A sizable Yukawa phase ip could generate values of the CP-asymmetry large enough to 
be relevant for future high-precision experimental analyses. However, a NNLO analysis of 
the theoretical prediction appears to be needed to reduce the presently large theoretical 
uncertainties and fully exploit such a measurement. 
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Figure 16: Constraints on |q| versus M H ± (in GeV) from B — > X s j, for <; u = 0.5 (left) and 
<; u = 1.5 (right). The white areas are excluded at 95% CL. In the upper panels, the phase has been 
varied in the range (p S [0, 2tt]. The lower panels assume real couplings. 
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Figure 17: Constraints on \<; u \ versus M H ± (in GeV) from B — > X s j, for ^ = 0. The white area 
is excluded at 95% CL. 



6. Discussion 

Imposing natural flavour conservation through discrete Z2 symmetries, one finds that the 
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Figure 18: Constraints on % versus <; u (95% CL) from B — > X s "f, assuming real couplings and 
taking M H ± = 150 GeV (left) and M H ± = 400 GeV (right). 



CKM phase is the only source of CP violation in the resulting 2HDMs. During the last 
thirty years, it has been common lore to assume that this is a more general fact, i.e. 
that the absence of tree-level FCNCs implies the absence of additional phases beyond the 
CKM one. The A2HDM provides an explicit counter-example, where FCNC couplings are 
absent at the Lagrangian level, while additional unconstrained complex phases generate 
new sources of CP violation. Since all Yukawa couplings are proportional to fermion 
masses, the A2HDM gives rise to an interesting hierarchy of FCNC effects, avoiding the 
stringent experimental constraints for light-quark systems and predicting at the same time 
interesting signals in heavy-quark transitions. The flavour-blind phases present in the 
model open a very interesting phenomenology which is worth to be investigated. The 
built-in flavour symmetries protect very efficiently the A2HDM from unwanted FCNC 
effects generated through quantum corrections. At the one-loop level the only allowed 
FCNC local structures are the two operators in ( 2.13| ), which could have very interesting 
(and computable) implications for B® mixing. 

Besides the fermion masses and mixings, the charged-scalar couplings of the A2HDM 
are fully characterized by three complex parameters In the previous sections, we have 
analyzed the impact of the contribution to different observables, where it is expected 
to be the dominant new-physics effect. Using conservatively estimated hadronic param- 
eters and up-to-date data, we have inferred the present constraints on the new-physics 
parameters involved in these processes. 

Leptonic tau decays provide a direct bound on the leptonic Yukawa coupling: |q \/M h ± < 
0.40 GeV" 1 (95% CL). From semileptonic processes constraints on the products <^\ U /M^ ± 
and ^*?rf/M^ ± are derived. The leptonic decays of heavy- light mesons allow us to disen- 
tangle the effects from q u and q. Thus, from B — > rv we derive an annular constraint 
in the complex plane ?*^/Mg- ± (figure [j]a), implying the absolute bound |?*?d/-M^ ± | < 
0.108 GeV~ 2 (95% CL). For real Yukawa couplings there is a two-fold sign ambiguity 
generating two possible solutions, the expected one around Ajj = (the SM ampli- 
tude dominates) and its mirror around = 2, corresponding to a new-physics con- 
tribution twice as large as the SM one and of opposite sign. The real solutions are 
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<^ d /Mf± G [-0.036,0.008] GeV" 2 and [0.065,0.108] GeV" 2 . 

Similar, but slightly weaker constraints on q*q u /M^j ± are obtained from the decays 
D — >■ fxu (figure [j]b) and D s — > (r, [i)v (figure ||); in the last case the bounds from B — ^ rv 
are used to get rid of the small ^ contamination proportional to the strange quark mass. 
The resulting absolute bound \tfs u /M^±\ < 0.6 GeV -2 (95% CL) is rather weak, but the 
upper limit corresponds to a new-physics contribution twice as large as the SM one, a very 
unlikely situation. The annular form of these constraints results in much stronger limits, 
once this possibility is excluded by other processes. For real Yukawa couplings, one finds 
sfSt/Affpt € [-0.005,0.037] GeV" 2 or [0.511,0.535] GeV~ 2 , at 95% CL. 

Owing to the quark-mass suppression, the absolute constraints obtained from lep- 
tonic decays of light mesons (figure ^) are obviously much weaker. However, the excellent 
experimental precision achieved in ir and K decays implies a narrow allowed annular re- 
gion. For real Yukawa couplings this translates into quite stringent bounds: q*^/M^r± G 
[-0.07,0.07] GeV" 2 or [8.14,8.28] GeV" 2 (95% CL). The uncertainties are dominated by 
the present theoretical knowledge of the ratio fx/ fir- 

Independent information is obtained from the semileptonic decays of pseudoscalar 
mesons, through the scalar form-factor contribution. One needs, however, to disentangle 
the dominant vector form-factor amplitude, which does not contain any charged-scalar 
effect and is correlated with the usual measurement of the corresponding CKM mixing 
factor. The present constraints from the ratio Hi{B — > Dtv t ) /Br(£? — > Dev e ), shown in 
figures | and |, are not very strong by themselves, but allow in combination with other 
processes the exclusion of the second real solutions in the Su^tf / M^ ± planes. A future 
measurement of the differential distribution in B — > Dtv t would obviously increase the 
sensitivity to the scalar contribution. In spite of the strange-mass suppression, the much 
higher experimental accuracy achieved in the analysis of K — > ttIv decays allows to derive 
the bound Re(^*? d /M| r± ) G [-0.16,0.30] GeV" 2 (95% CL). This already excludes the 
second real solution (a scalar amplitude larger than the SM one) obtained from K^l^^i- 

Combining the information from all leptonic and semileptonic decays analyzed, one 
gets the constraints shown in figure || 

The flavour-conserving decay Z — )• bb provides a very stringent constraint on \q u \. Since 
Vtb ~ 1, the one- loop contributions involving virtual top quarks completely dominate both 
the SM (W 1 * 1 ) and the new-physics (H^) radiative corrections. In contrast to leptonic 
and semileptonic processes, where the charged-scalar effects are necessarily proportional 
to q, the Z — » bb amplitude gives direct access to ? u and <id- Owing to the relative 
factor rrib/mt which suppresses the (>d contribution, one gets finally the constraints on 

shown in figure |9] (assuming |q| < 50). At 95% CL, we obtain \q u \ < 0.91 (1.91), for 
M H ± = 80 (500) GeV. The upper bound increases linearly with M H ± , implying \ q u \/M H ± < 

0.0024 GeV" 1 + < 0.011 GeV" 1 , where we have used the LEP lower bound on the 

M H ± ' 

charged-scalar mass M H ± > 78.6 GeV (95% CL) [jl], |8l] ]. Together with the tau-decay 
constraint on |q|/M#±, this gives the limit |? U ^*|/M^ ± < 0.005 GeV -2 , which is much 
stronger than the information extracted from the global fit to leptonic and semileptonic 
decays. 

Quite similar information can be extracted from B° mixing, which is also dominated 
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by one-loop contributions involving virtual top quarks. The smallness of the m s /M\y 
ratio implies that SU(3)-breaking corrections are negligible; therefore, the charged-scalar 
contributions cancel in the ratio Am B o / Am B o , which can be used in the CKM fit. Only two 

s d 

AB = 2 four-quark operators are numerically relevant; the one generating the leading SM 
amplitude gets new-physics contributions proportional to |?«| 2 ' 4 , while the other operator 
generates subleading corrections proportional to (<j*<j d ) 1,2 m^/M^ / . Scanning the parameter 
ranges |q| < 50 and (p £ [0,2tt], where <p is the relative phase between s u and q, the 



measured B® mixing amplitude implies the constraints shown in figure 10, in the plane 
M H ±-\<; U \. At 95% CL, one gets | ? u | < 0.00279 M H ± +0.27+117/A%± , for M H ± G [80,500], 
in GeV units. 

The charged-scalar contribution could accommodate a large B® mixing phase, without 
spoiling the agreement in the B^ system, although it is not possible to reach a value as 
large as hinted at by the present DO central value, which is however at odds with the rate 
difference being unaffected by new physics (unless the calculation of the rate difference 
is affected by problems regarding the OPE). If confirmed, a large phase <j> s would point 
towards large values of |^|, small charged-scalar masses and a sizable Yukawa phase (p. 
The preferred negative sign for the asymmetry would require ip E [tt/2,tt], [3tt/2,2tt]. 
Additional contributions to (j) s could be induced by neutral scalar exchanges, through the 



effective FCNC operator in Eq.(2.13) appearing at the one-loop level. Large Yukawa phases 
could be constrained by other CP-violating observables not yet included in our analysis. 
A detailed discussion of these effects is postponed to future work. 

The observable ex leads again to a similar constraint, even slightly more restrictive 
than the ones from B° mixing and Z — > bb. Although CP violating, this observable is 
insensitive to the new-physics phases, as the relevant contributions involve \s u \, only. We 
obtain at 95% CL < 0.560 + 2.647 W~ 3 M H ± - 1.049 lO^M 2 ^ + 6.153 l(r 10 Mj^ ± in 
units of GeV. 

The radiative decay B — > X s ^f provides another important source of information. There 
are two different charged-scalar contributions, proportional again to \q u \ 2 and but in 
this case the two have similar sizes. Their combined effect can be quite different depending 
on the value of the relative phase ip. This results in rather weak limits because a destructive 
interference can be adjusted through this phase. The resulting constraints on and |<^| 



are shown in figure 15, varying the charged-scalar mass in the range M H ± G [80, 500] GeV. 
Scanning the phase (p in the whole range from to 2ir, and imposing < 3, one finds 
roughly < 20 (95% CL). Much stronger bounds are obtained at fixed values of the 

relative phase. Assuming real values of s u and q (i.e. ip = or 7r), one finds that couplings 
of different sign are excluded, except at very small values, while a broad region of large 
equal-sign couplings is allowed, reflecting again the possibility of a destructive interference. 



Figures 16, 17 and [L8| show the sensitivity of the B — > X s j constraints to the different 
unknown parameters: M H ±, \q u \, |^| and (p. 

The constraints discussed so far apply to the general A2HDM framework, with three 
arbitrary complex parameters The limits become of course much stronger in particular 



models where these parameters are correlated. Figures 19 show the combined constraints 
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Figure 19: Constraints on Mjj± (in GeV) versus tan/3 (95% CL), in the 2HDM models of types 
I (upper-left), II (upper-right), X (lower-left) and Y (lower-right). 



on the tan f3~M H ± plane for the different Z2 models. The bounds from Z — > bb, ck, Am#o 
and B — > X s ^f are obviously identical for the models of type I and X and also for type II 
and Y. In the type I/X case, ^ = ^ = q u <^d = cot 2 (5 and the scalar amplitudes grow for 
decreasing values of tan/3. For type II/Y, this behaviour is only observed in the ^ term, 
while ?| = tan 2 /? and s u Sd = — 1; the decay B — > X s j provides then a very strong lower 
bound on the scalar mass, M H ± > 277 GeV (95% CL), due to the constructive interference 
of the two contributing amplitudes. The q coupling gives rise to different constraints from 
leptonic and semileptonic decays in each of the four models. Our results agree with the 
qualitative behaviour found in previous analyses Jl 



38, 40, 41 



136, [L45L 146, 147, 



148] , 149| , |150| , the small differences originating from the slightly different inputs adopted. 
The A2HDM is not the most general version of a 2HDM without treedevel FC- 
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NCs. To avoid the unwanted FCNCs one just needs diagonal Yukawa matrices Yt in 
the fermion mass-eigenstate basis, i.e. Y d = diag(y rf , y s , y b ), Y u = diag(y u , y c , y t ) and 
Yi = diag(y e , ?/«, y T ), with arbitrary parameters y{. This more general scenario can be 
formally described by the Lagrangian ( 2.1 0| ) with the substitution ?/Mj — > Yf. One could 



still use nine dimensionless parameters <jj = yj/rrif, one for each charged fermion [151], but 



in this case this is just a redefinition of the Yukawa couplings yj because a priori nothing 



relates them to the fermion masses [152]. The hierarchy of couplings characteristic of the 
A2HDM ansatz is lost and one can no longer justify that the leading charged-scalar effects 
originate in the heavier fermion couplings (it becomes an assumption). With this caveat 
in mind, our results can still be applied in this case, but most correlations among different 
processes disappear because the associated constraints correspond now to different $f pa- 
rameters. For instance, the constraint in ([0|) refers to y^r^f and figures [l] to s*% (left) 
and ?*? c (right). 

The A2HDM provides a general setting to discuss the phenomenology of 2HDMs, 
satisfying in a natural way the requirement of very suppressed FCNC effects. The align- 
ment conditions imply Yukawa couplings proportional to the corresponding fermion masses, 
which is supported by the data (bounds of order 1 for the parameters). While including 
as limiting cases all Z2 models, the A2HDM incorporates possible new sources of CP vi- 
olation through the phases. The additional freedom introduced by these phases makes 
easier to avoid some low-energy constraints, resulting in weaker limits than in the usual 
scenarios with discrete Z2 symmetries. A detailed analysis of CP-violating observables 
is clearly needed to investigate the allowed ranges for these phases and their potential 
phenomenological relevance [35]. 



At the moment, the data does not show any clear deviation from the SM. Therefore, 
we have derived upper limits on the Yukawa parameters. Nevertheless, we have already 
pointed out that the A2HDM could account for a sizeable B® mixing phase, as suggested 
by the present B s — > J/tjxf) and like-sign dimuon data. Our bounds could be made stronger, 
adopting more aggressive estimates for the hadronic parameters entering the analysis, but 
we have preferred to be on the conservative side and infer solid limits for later use. Im- 
provements are to be expected on one hand from better theoretical determinations of the 
hadronic inputs, and on the other hand from more accurate measurements at NA62 (kaons), 
LHCb (Am^ )S , B s — > J/ificf)), a future Super-B factory (t, b — > 57, Am^, B — > tu, B — >• Dlv\ 
or a linear collider with Giga-Z option (Rb). The agreement of the different bounds in the 
vicinity of zero is trivial, when the SM agrees with the data. If signals for new-physics are 
found at LHC, the analysis presented here will be capable of quantifying the agreement 
(or disagreement) of the data with the A2HDM, and with the different implementations of 
the 2HDM based on Z2 symmetries, in one step. 

Note added 

After this work was finished, two relevant papers have been posted in the archives. In 
Ref. fll53j [ an approximate solution to the renormalization-group equations of the A2HDM 
is analyzed and the generated FCNC terms are studied numerically; the results presented 
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there agree with our FCNC operator ( |2.13| ) and it is concluded that the induced FCNC 
effects are well below the present experimental bounds. Ref. [154] analyzes the strength 
of FCNC effects mediated by neutral scalars in a minimal-flavour-violating framework 
containing two Higgs doublets, assuming a perturbative expansion around the type II 
model. The tree-level alignment conditions of Ref. |TJ are reproduced, the one-loop FCNC 
structures in ( p. 13 ) are discussed and their coefficients are estimated at large tan (3 in the 
decoupling limit. The phenomenological analysis of Ref. [154] emphasizes the potential 
relevance of the flavour-blind phases present in the A2HDM to accommodate the recent 
hints of a large B® mixing phase through neutral-Higgs exchange. 
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A. AF = 2 effective Hamiltonian 



A.l AB = 2 



At lowest order, the AF = 2 transitions are mediated by box diagrams with exchanges of 
and/or H^ 1 propagators. Performing the matching between the A2HDM amplitude 
and the low-energy effective Hamiltonian T-L^f =2 , at the scale fitw ~ M\y,m t , one obtains 
the Wilson coefficients Ci(/i). We have derived the LO results given in table where 
xy/ = ml/M^Y and xu = m t/Mj I± . They can be expressed in terms of the two four-point 
functions J155| : 



D (m 1 ,m 2 ,M 1 ,M 2 



m. 



I log 



(ml 



+ 



ml){m% - Mf) (ml - Mf) 



+ 



Ml log 



Ml 



(Mf - m{)(M( - ml)(M( - Mf ) 



Ml 



log 



+ 



Mf 



(Ml-mj)(Ml 



ml)(Ml-Ml) 



(A.l) 
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4 



- log (5) 

Lhimi - "' 2 " U ' ■ U2i " { m \-rn\){rn\-MUrn\-Ml) + 

f4 l„„ 



+ (Mf - m\){Ml - m\)(Ml - Ml) + 



M 2 4 log 1^ 



(M$-ml)(M$-ml)(M$-M?) 



(A.2) 



through (i = 0, 2) 



A(m,Mi,M 2 )= lim A(m,m 2 ,Mi,M 2 ) , (A.3) 

m2->m 

DAm,M)= lim DJm, M, M 2 ) , (A.4) 

M 2 ->M 

D 2 (m,Mi,M 2 ) = Z} 2 (m,Mi,M 2 )-Z} 2 (0,Mi,M 2 ). (A.5) 

These one- loop contributions involve virtual propagators of up-type quarks (u,c,t). 
Once the GIM cancellation is taken into account, the up and charm contributions vanish 
in the limit m U)C — > 0, which we have adopted. Thus, the B meson mixing is completely 
dominated by the top-quark contributions (the different CKM factors have all a similar 
size for B° d mixing, V* d V ub ~ V* d V cb ~ ~ AA 3 , while in the B° case V* s V ub ~ ,4A 4 

and V^V^ ~ ^s^, ~ AA 2 ). Since the scalar couplings are proportional to quark masses, 
we have maintained the masses of the external down-type quarks. In the limit — > 0, we 
reproduce the results given in [121]. The only Wilson coefficients which are not suppressed 
by powers of are Cvll and Cg RR . Therefore, for all practical purposes, one can neglect 
the remaining operators. 

The running for Of RR is performed using the results of ] |119[| , 

(A.6) 



CsrrO^) I = | b7ll(Mb)]sRR [?7120&)]SRR j ( Cg RR (iw) 
C lKK^b) ) \ [f?2l(W>)]sRR [??22(W>)]SRR / V CsKR&tw) 



with 

failO"6)]sRR = 1.0153 r?" - 6315 - 0.0153 rj° 5 71M , (A.7) 

[vM]srr = 1-9325 (r^ - 6315 - r,° 5 7184 ) , (A.8) 

feiOMsRR = 0.0081 (4 718i - % - - 6315 ) , (A.9) 

[r?22(^)] S RR = 1.0153 % a7184 - 0.0153 r^ 06315 . (A.10) 

These are leading-order expressions, but they have been evaluated with the two-loop ex- 

pression for a s in 775 = 013 (5 ) ~ 0.7. 

«s («>) 

The corresponding matrix elements are given by 

(0 VLL )(M)= ~m B o/| o B VLL 0u), (A.ll) 

3 d d 



(Oi srr )(m) = ~ (—rrr-n ) "bs&sPV) > ( A - 12 ) 

24 \m b {n) + m d (n) 1 



(02 SRR )(M) = -5 ( , ™f p r V m ^^°g 2 SRR (M) , (A-13) 
2 \m b (n) +m d {fi)J d *d 



- 34 - 



Oi 

qVRR 

o 2 LR 



SLL 
SLL 



qSRR 
Q SRR 



Ci{ntw) 

(4x w + x 2 w )M^D 2 (m t , M w ) - 8x 2 w M^D (m t , M w )+ 

+ 2\s u \ 2 x 2 w [M 2 ,D 2 (m t ,M w ,M H ±) - 4M^D (m t ,M w ,M H ±)] + 

+ \c; u \ 4 x 2 w M 2 v D 2 (m t ,M H ±) 

~\q d \ A x H M 2 v D 2 (m u M H± ) + \q d \ 2 M^D 2 (mt,M w ,M H ±)] 
xw [\c;d\ 2 K\ 2 M^D 2 {m t ,M H ±) + 2Re(^u)M^D 2 (m u M w ,M H ±)] 
[^d\ 2 K\ 2 xwM^D {m t ,M H ±)-A\^ d \ 2 M^D 2 (m t ,M w ,M H ±)+ 
+ (M 2 + \<, u \ 2 )xwM^D 2 (mt,M w , M H ±)] 

2 [(^* d ) 2 M^D (m t , M H± ) + 2^*M^D (mt, M w , M H± )] 



r, m d rn b 

r, rn d in b 
w 



[(^*) 2 M^D (m t ,M H± ) + 2^*M^D (m t ,M w ,M H ±)] 



Table 4: Leading-order Wilson coefficients for the AB = 2 operators given above. The quark 
masses from the scalar couplings are to be taken at the matching scale fj-tw- 



the Bi(n) parametrizing the deviation from the naive factorization limit. These Bi(fi) 



factors have been evaluated in the quenched approximation on the lattice in [156], using 



a different operator basis. The connection reads (see again [11£], given here with both 
operators in the same scheme) 



(A.14) 



From [156 1 we arrive at the values given in table || by adding again all systematic uncer- 
tainties linearly. 







B° 

^ s 


B^ IS (m b ) 
BfS(m b ) 


0.83 ±0.03 ±0.06 
0.90 ±0.06 ±0.12 


0.84 ±0.02 ±0.06 
0.91 ±0.03 ±0.12 



Table 5: £?-parameters for B d mixing from [156|. Systematic errors added linearly. 



The wanted B° d -B Q d mixing amplitude is given by 



0|-ivAB=2|^0\ 



(B»\H 



G 2 M 2 



off 



"2 - 

-B B o i]b(xw,xh) C V ll(m*wO+ 
6 d 



(A.15) 
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(m 6 (/x 6 ) + m d (n h )y 
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with 

/ _ 5 

Bsrr(m&) = 



T2 B ^) ] (A.16) 



§-E>3,d(^f>) — §-E>2,d(w> 



From this, we get the relevant observables as 



Am B° = -^\(B° d \U* B = 2 \B° d )\ , (A.17) 



d m R o 



V = -Arg [(B d \Ht e B = 2 \B° d )] . (A.18) 



The analogous expressions for B®-B® mixing are trivially obtained changing the label 
d to s everywhere. 

A.2 AS = 2 

For the Kaon mixing amplitude, we have calculated the LO matching coefficients completely 
anologous to the AB = 2 coefficients, keeping the charm mass finite. Due to the strong 
suppression of all other operators by light quark masses we can choose the LO matching 
coefficients to be 

C£ VLL = (4x c ^ + x c ^ 2 )M^D 2 (m c , M w ) - 8x c ^ 2 M^D (m c , M w ) , 

CcVll = ^ x w + x c ^ 2 )M 2 v D 2 (m c , m t , M w ) - 8x c ^ 2 M^D Q (m c , m t , M w ) + 

+2| ?u | 2 x^ 2 [M^D 2 {m c ,m t ,M w ,M H ±) - 4M^D (m c ,m t , M w , M H ±)] + 

+\q u \ A x c ^ 2 M 2 v D 2 {m c ,m t ,M H± ) , (A.19) 
C Ovll = ( 4x w + x 2 w )M^D 2 (m t ,M w ) - 8x 2 w M^D {m t ,M w ) + 

+2\<; u \ 2 x 2 v [M^D 2 (m t , M w , M H ±) - 4M^D (m t , M w , M H ±)] + 

+ \q u \ i x 2 w M^D 2 {m t ,M H ±) , 
C Qi =0 (i + VLL) , 

where the loop functions Dq j2 have been defined in appendix A, and = m c m t /M'^ v . In 
the calculation, we use the NLO results for the SM which have been calculated in |157| , |158| . 
while the NLO charged scalar contributions to the top contribution are again taken from 
1 121], corrected and applied to our scenario. 
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